Landau Damping of Baryon Structure Formation in the Post Reionization Epoch 
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It has been suggested by Chen and Lai that the proper description of the large scale structure 
formation of the universe in the post-reionization era, which is conventionally characterized via gas 
hydrodynamics, should include the plasma collective effects in the formulation. Specifically, it is 
the combined pressure from the baryon thermal motions and the residual long-range electrostatic 
potentials resulted from the imperfect Debye shielding, that fights against the gravitational collapse. 
As a result, at small-scales the baryons would oscillate at the ion-acoustic, instead of the conventional 
neutral acoustic, frequency. In this paper we extend and improve the Chen-Lai formulation with 
the attention to the Landau damping of the ion-acoustic oscillations. Since Te « Ti in the post- 
reionization era, the ion acoustic oscillations would inevitably suffer the Landau damping which 
severely suppresses the baryon density spectrum in the regimes of intermediate and high wavenumber 
k. To describe this Landau-damping phenomenon more appropriately, we find it necessary to modify 
the filtering wavenumber fc/ in our analysis. It would be interesting if our predicted Landau damping 
of the ion-acoustic oscillations can be observed at high redshifts. 

PACS numbers: 



I. INTRODUCTION 



The evolution of the cosmic structure has been an im- 
portant subject in cosmology. In the post-decoupling era 
the neutral atoms (the baryons) ceased to feel the pres- 
sure from the photons. So the only force that fought 
against the gravitational pull was its own thermal pres- 
sure. Due to the Hubble expansion, the baryon temper- 
ature gradually decreased and so was its thermal pres- 
sure. The situation changed drastically at the subsequent 
reionization epoch, during which baryons turned into 
fully ionized plasmas and were heated up to T ^ lO^'i^T by 
the photo-ionization process Such reheating process 
raised the Jeans scale to the size of the galaxies [2| and as 
a result the course of the structure formation was signifi- 
cantly altered. Traditionally the study of the baryon evo- 
lution has been conducted by staying within the gas hy- 
drodynamics framework while changing the mean molec- 
ular weight between the pre- and the post-reionization 
epochs to account for the effects of ionization. However, 
it is known that gas hydrodynamics, which assumes local 
charge neutrality, cannot exhibit crucial collective phe- 
nomena that is unique to plasma, such as Landau damp- 
ing. It is therefore natural to wonder whether relevant 
plasma collective effects in cosmic evolution might have 
been overlooked. 

In 2007, one of us (Chen) and Lai first considered the 
plasma effects on the baryon structure formation in the 
post-reionization epoch, where the conventional baryon 
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acoustic waves at small scales are replaced by the ion 
acoustic waves in plasmas [4|. A new formulation based 
on the Maxwell-Einstcin-Boltzmann equations for pure 
hydrogen ions was introduced in that paper. To obtain 
the baryon power spectrum, the authors employed the in- 
stantaneous Jeans wavenumber kj and interpolated the 
solutions from the opposite limits of the wavenumber k. 
Interesting as it is, this new plasma hydrodynamic for- 
mulation has some shortcomings, however. In particu- 
lar, the spectrum in the intermediate regime of k, which 
is the most crucial regime of interest, was unresolved. 
To properly account for the impact of the Hubble ex- 
pansion, it was suggested that the relevant scale should 
be the integrated and smoothed 'filtering wavenumber' 
kf rather than the instantaneous Jeans wavenumber kj 
[a, IB]- Motivated by the reasons, in this paper we analyti- 
cally re-solve the baryon density spectrum under the Lan- 
dau damping from the corrected plasma hydrodynamic 
equation of hydrogen-helium plasmas. Having the solved 
spectrum, we are able to investigate the impacts of the 
Landau damping effect on the evolution of baryon over- 
density and on the filtering wavenumber kf due to the 
decaying amplitude of baryon oscillation modes. 

The layout of the paper is the following: In SecHIl 
we briefly review the concept of ion acoustic wave and 
its associated new sound speed for baryons. In SeclIIII 
we introduce the Landau damping of the baryon density 
spectrum. This is carried out through the analytical so- 
lution of the plasma hydrodynamics equation. Then in 
Sec lIVI we investigate the influence of Landau damping 
on the baryon spectrum and derive a modified filtering 
wavenumber fc^ in place of fc/. Finally we compare the 
spectra obtained from plasma and gas approaches, re- 
spectively, and discuss the possible ways to distinguish 
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the two models. respectively, for the neutral gas. In general, the i on flu id 

leads to a lower sound speed by a factor of \/ lij [i ^ 
1.12, which corresponds to ^ 10% increase of the Jeans 
II. ION ACOUSTIC FORMULATION wavenumber fcj. 



An ordinary sound wave in a neutral gas is a pres- 
sure wave that propagates through the thermal collisions 
among neighboring particles. Without collisions, the 
sound wave would not occur. However for an ion acoustic 
wave in a fully ionized plasma, an extra repulsive, non- 
collisional, electrostatic pressure exists due to the im- 
perfect Debye shielding at finite temperature. Hence an 
ion acoustic wave can travel through the medium without 
collisions 0, Q . Based on the Chen-Lai formulati on [1 , Q 
and the fluid equations for matter overdensites (5l.ll0l|. we 
obtain the coupled fluid equations for hydrogen-helium 
baryon (ignoring the electron mass): 

liu + 2Hhik = 4^Gp {hhk + /d<5dfc) , (la) 
'k + 2i/(5b + - 4^Gp (/fcJfcfe + /d(5rffc) , (lb) 

where the subscripts h and d indicate the baryon and the 
dark matter, respectively, (5 = ipj p is the fractional over- 
density, /(, = f2b/f2m and = 17c; /fim arc the baryon 
and dark matter fractions and Cg is the sound speed for 
ionic baryons. According to plasma physics, the sound 
speed is defined by 



cl = dPb/dpb+S'i>ern/5pb, 



where 



Spb = mnSnH + rriHeSnHe, (2a) 
SPb = JbkBT{6nH + SuHe), (2b) 

<5$em = JekBT{6nH + ZSuHe), (2c) 

Ze is the helium ion charge and $e„i is the combined 
residual electrostatic potential of hydrogen and helium. 
Here we have assumed the temperature Te ~ Th ~ 
Thc ~ T because the time scale for H-He plasma to 
reach thermal equilibrium is short compared with that 
of the duration of the reionization epoch [ll| . In a typ- 
ical plasma 7e = 1 for isothermal compression condition 
is often assumed since the thermal velocities of electrons 
were ~ ^fnh^Jm'e times larger than the fluid sound speed, 
which enables electrons to interact with each other over 
several wavelengths within one oscillation Q . Hence in- 
voking the primordial helium abundance Y ~ 0.24 [l2j 
deduced from the big bang nucleosynthesis (BBN), we 
can express the ion sound speed Cg in terms of the ef- 
fective mean molecular weight /i, which combines 7e and 

76, 



3 \irrip ' 



(3) 



with p. = 0.75 for Z = 2 and p. = 0.77 for Z = 1. Note 
that in contrast p — 0.59 and 0.62 for Z = 2 and 1, 



III. BARYON DENSITY SPECTRUM WITH 
LANDAU DAMPING 

It is known that plasma waves suffer the Landau damp- 
ing when the thermal velocity of the plasma particles is 
slightly below the phase velocity velocity. Residing on 
the down-slope of the Maxwellian particle velocity dis- 
tribution, the plasma wave tends to lose energy to more 
particles and gain energy from less particles. This results 
in the decline of the plasma wave amplitude. Learning 
from that, we expect the similar effect on the small-scale 
baryon oscillations, the ion acoustic modes. Mathemati- 
cally, Landau damping is a consequence of the singularity 
w — k • V in the Vlasov equation, which is absent in the 
fluid theory. To account for the Landau damping effect 
in the fluid model, one may separately calculate the Lan- 
dau damping rate and then multiply it to the imaginary 
part of the plasma fluid solution. 

To analytically study the baryon density spectrum un- 
der the Landau damping effect, we start with the coupled 
fluid equations, Eq. ([Ta|) and Eq. (jlb[) . with fb — and 
fd ~ ^ for simplicity. Replacing the variable t by a via 
the relation t cx a^^^ in the matter-dominant era, we ar- 
rive at 
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(4a) 
(4b) 



where kj is the co-moving Jeans wavenumber of the ionic 
baryon. The solution of Eq. (|4ap gives Sdk oc D-^-(a) cx a 
when dropping the decaying mode, while the solution of 
Eq. (j4bp depends on the thermal history of reionization. 
We adopt the temperature evolution of the intergalactic 
medium (IGM), T oc a"°-^®, which rcsuks from the bal- 
ance between the photo-heating process and the Hubble 
expansion cooling process [l3|. In fact the recent mea- 
surements indicate a second rehea ting process at z ~ 3 
due to the ionization of Hen [13, Il5| . To simplify the 
case, we assume a single reheating process and the tem- 
perature evolution T = Tqo"^ which leads a constant 
Jeans wavenumber fcj(a) = fcjion, where the subscript 
ion stands for the ending time of the reionization pro- 
cess. The equation with arbitrary thermal evolution can 
be analytically solved via Green's function method jlGj . 

Inserting the solution of Eq. Pa|) . Eq. (j4b[) can be 
rewritten as 
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(5) 



A = S/Sion and k = k/kjion- Note 
a. Since we are interested in the relation 
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between the baryon and the dark matter evolutions, it 
is convenient to express the baryon density in terms of 
that of the dark matter. Using the initial conditions: 1) 
Afcfeion = ^dkion = 1 and 2) dAbk/da = dAdk/da = 1 at 
diom we obtain the normalized baryon density spectrum: 
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(6) 



where Q = vT^24k^ and 
5 + Q k2 



Ai = 



2Q 1 



A. = 



2Q 1 



As > 24, Q becomes imaginary. We may therefore 
rewrite Eq. (j6|) as 



A 



bk 
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5/4 



ia,K^ > 1/24) 
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in f — In a I +cosf — lna| 
V4 J V4 / 



bi{a, k) + 
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(7) 



where q = \/2Ak^ — 1 and bi(a, k) denotes the oscillatory 
part of Abk/Adk- If k — > oo or equivalently q — > oo, 
Abk/Adk asymptotically approaches a~^/^ cos(g/4). Wc 
see the amplitude of the baryon oscillations decreases as 

Because Te « Ti, the thermal velocity of the ions is 
close to the sound speed. We then expect that the Lan- 
dau damping would occur to the baryon fluid oscillations. 
By invoking the nonlinear empirical formula which re- 
lates the real and the imaginary parts of the ion acoustic 
oscillation frequency , we introduce the Landau damp- 
ing effect to Eq. ([7|) by multiplying the oscillation term 
with an exponential factor, 



^(a,^2 > ^/24) = e-(«'-'/4"'-)'""6,(a,^) + 
Adk 1 -I- 

where 
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0.4. 



(9) 



This means the ion acoustic wave would be Landau- 
damped to one e-fold of its initial amplitude within just 
~ 0.4 of one oscillation period. Such severe damping ef- 
fect causes the baryon oscillations to be totally smoothed 
out at late times. 



IV. RESULTS 

With the Landau damping introduced, we now calcu- 
late the complete baryon density spectrum by combining 
Eq. ^ and ([5]). Figure [T] shows the result (in solid curve) 
at z = 3 with Zion = 7. The dashed and dotted curves 
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FIG. 1: The baryon density spectrum under Landau damp- 
ing(solid) at z = 3 with Zion = 7, superimposed by the dashed 
and dotted curves standing for the results from gas liydro- 
dynamics (/i) and that from plasma physics before applying 
landau damping (jl) respectively. 
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FIG. 2: The same plot as Fig. [T] but has the dotted curve 
represent the fitting function Abk/^dk = exp(— fc^//cj^) with 
the modified filtering wavenumber k'f. 



correspond to the results obtained from the neutral gas 
hydrodynamics and plasma physics without applying the 
Landau damping. We see that the dotted curve shifts 
from the dashed one by ~ 10% due to the increase of the 
Jeans wavelength mentioned in Sec. |TT1 When suffering 
the Landau damping, the modes of baryon overdensity at 
large k diminish, leading to a suppression of the baryon 
density spectrum in the intermediate regime of k. Such a 
suppression, however, accidentally compensates the shift 
of the spectrum. Thus the filtering wavenumber kf pro- 
posed by Gnedin and Hui Q to characterize the baryon 
structure would not be suitable when the Landau damp- 
ing effect is included. 

In an expanding universe, the Jeans scale changes with 
temperature and therefore with time. On the other hand, 
it takes about one Hubble time for baryons to respond to 
the changed pressure. Therefore the filtering wavenum- 
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FIG. 3: The comparison of two characteristic scales in gas 
and plasma approaches versus redshift z. The solid curve is 
k'j /k,j and the dashed curve is kf /kj where kj is the ordinary 
Jeans scale for baryon gas. 
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FIG. 4: The comparison of baryon density spectra at z — 6.5 
and z = 5 for gas (dashed) and plasma (solid) approaches. 



ber kf was introduced @, @ to account for the migration 
of the Jeans wavcnumber over the time scale relevant to 
the acoustic oscillations, which therefore depends on the 
entire previous thermal history of the universe j[l7l . In 
the matter-dominant era fc/ is related to kj by 



kj 



da' 



1 - 



1/2- 



(10) 



When the baryon fluctuations suffer the Landau damp- 
ing effect, a correction to fc/ must be taken. The original 
kf was deduced by supposing the solution Atk/^dk = 
1 — k'^/kj sa exp{—k'^/kf) under the small k approxi- 
mation (fc — > 0). With the exponential damping fac- 
tor included, we find it reasonable to directly assume 
Ahk/ Adk{ct, k) « exp{—k'^/kj), where kf is the modified 



filtering wavcnumber. Assuming kf /kj = cr, we have 

exp(-4) = e-°-4('/4)'""5,(a,«) + (11) 

Upon taking logarithm and then square-root on both 
sides, we obtain 



^7 



The right-hand side of the above expression approaches 
a constant for k G [0.8, 2] at z < 5, which falls in the 
intermediate regime of k. Hence the modified filtering 
wavcnumber k / can be obtained by substituting an arbi- 
trary K within this range. Figure [2] is the same plot as 
FigUbut with an additional dotted curve that illustrates 
the well-fitted function of exp(— fc^/fcj). In Fig. [3] we 
compare the two characteristic scales kf (dashed) and kg 
(solid) at different redshifts. The two curves converge to 
each other at low redshifts, while in general kf is smaller 
than kf until 2^1. 



V. CONCLUSION 

In summary, the baryon overdensities has three stages 
of evolution: 1) before the photon-baryon decoupling, 
where baryons were in a plasma state but were tightly 
coupled to and dominated by the photons; 2) during the 
dark age, where baryons formed neutral gases and grew 
with dark matter overdensities; and 3) after reionization, 
where baryons turned into a plasma again and followed 
the plasma hydrodynamics. Although before the decou- 
pling baryons were in a plasma state, the photon-baryon 
fluid was dominated by the photon pressure rather than 
the combination of the baryon thermal and electrostatic 
pressures. This resulted in a sound speed that is much 
larger than the ion thermal velocity. The plasma ef- 
fects were negligible in that epoch. This aspect has been 
widely recognized. 

What appears less recognized is the influence of the 
plasma collective effects in the post-reionization epoch. 
Considering the plasma effects during this epoch, we 
point out that the baryons would experience two major 
impacts: the change of the sound speed and the Landau 
damping of the baryon oscillations. Numerically speak- 
ing, the first effect is relatively minor compared with 
the second. Since Ti ss T^, we expect that the Landau 
damping would be quite effective in dissipating the am- 
plitude of baryon oscillation modes. This might influence 
the subsequent galaxy formations history. According to 
Eq. ([7]), the amplitude of baryon oscillations scales as 
Q,-5/4^ Therefore the discrepancy between the gas and 
the plasma approaches is most significant in the early 
stage of structure formation immediately after the reion- 
ization. Figure 2] is a comparison of the baryon density 
spectra based on the gas (in dashed line) and the plasma 
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(in solid line) approaches at different redshifts {z = 6.5 
and z = 5) with Zio„ = 7. If the baryons behave as a 
neutral gas, then we expect a sizable signature of baryon 
acoustic oscillations in the matter power spectrum below 
the scale k > 0(l)/iMpc~^ after the photo-reheating pro- 
cess. Conversely, if the baryons act more like a plasma, 
then wc should find little sign of oscillations but instead 
the Landau damping of it, whose rate maybe measured 
via tracing the spectrum back to the early time. So far 
the observations of the matter power spectrum have not 
yet reached the sensitivity for such small scales at high 
redshifts [l^. It would be interesting if this regime can 



be explored in the near future, so that we can learn more 
about the nature of the cosmic structure formation. 
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